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ABSTRACT 

The stability of a rotating fluid disk to the formation of spiral arms is 
studied in the tightwinding approximation in the linear regime. The dispersion 
relation for spirals that was derived by Bertin et al. is shown to contain a 
c$ ! new, acoustic instability beyond the Lindblad resonances that depends only 

on pressure and rotation. In this regime, pressure and gravity exchange roles 
as drivers and inhibitors of spiral wave structures. Other instabilities that 
are enhanced by pressure are also found in the general dispersion relation by 
including higher order terms in the small parameter 1/kr for wavenumber k 
and radius r. We identify two important dimensionless physical parameters: 
e = 2nGao/(rK 2 ), which is essentially the ratio of disk mass to total mass (disk 
and halo), and a/(nr), which is the ratio of epicyclic radius to disk radius (<t 
is the mass column density, k is the epicyclic frequency, and a is the sound 
speed). The small term ( = (k 2 r 2 + m 2 ) is an additional parameter that is 
purely geometrical for number of arms m. When these terms are included in the 
dispersion relation, the oscillation frequency becomes complex, leading to the 
growth of perturbations even for large values of Toomre's parameter Q. The 
growth rate is proportional to a linear combination of terms that depend on e 
and a/(nr). Instabilities that arise from e are termed gravitational- curvature 
instabilities because e depends on the disk mass and is largest when the radius 
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is small, i.e., when the orbital curvature is large. Instabilities that arise from 
a/(nr) are termed acoustic- curvature instabilities, because they arise from only 
the pressure terms at small r. 

Unstable growth rates are determined for these instabilities in four cases: 
a self-gravitating disk with a flat rotation curve, a self-gravitating disk with 
solid body rotation, a non- self-gravitating disk with solid body rotation, and 
a non-self-gravitating disk with Keplerian rotation. The most important 
application appears to be as a source of spiral structure, possibly leading to 
accretion in non-self-gravitating disks, such as some galactic nuclear disks, disks 
around black holes, and proto-planetary disks. All of these examples have short 
orbital times so the unstable growth time can be small, even when only terms 
of order e contribute. 



1. Introduction 

Spiral galaxies are characterized by bright "arms" spiraling out from a region near the 
center. Differential rotation will shear and wind these arms quickly if they are material 
features, so Lindblad (1958) and Lin & Shu (1964) developed a theory of density waves to 
overcome this winding dilemma. Lin & Shu (1964,1966) also obtained the dispersion relation 
for these waves, which is the relation between frequency and wavenumber. An important 
discriminant in this dispersion relation is the stability parameter Q for axisymmetric 
disturbances (Toomre 1964); when Q > 1, the disk is stable against ring-like disturbances. 

Lau & Bertin (1978) included additional terms that treated tangential forces for fluid 
spiral waves in a uniform disk, finding an additional destabilizing term they called J. 
They used a WKB approximation and ignored curvature terms, which scale inversely with 
galactocentric radius. Goldreich & Lynden-Bell (1965), Zang (1976), and Toomre (1981) 
also studied azimuthal forces, by considering the temporal response of shearing wavelets. 
Toomre (1981) termed the mechanism responsible for the spectacular growth of shearing 
waves a " swing amplifier" . He found that spiral waves can grow for a short time even when 
Q > 1, as long as Q is not too large. 

In section |3.2| below, we discuss a new instability in the usual spiral wave equations 
derived by Bertin et al. (1989; hereafter BLLT) that is relevant beyond the Lindblad 
resonances, i.e., inside the ILR and outside the OLR, even when Q > 1. This is a regime 
that BLLT did not consider. The new instability depends on shear and self-gravity as in 
the BLLT derivation, but it also has a component in the absence of self-gravity that arises 
only from pressure and rotation. We therefore refer to it as an acoustic instability. 
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We also derive dispersion relations for fluid disks considering the curvature terms and 
other terms that were ignored in these previous studies, such as radial variations of the 
basic properties of the disk. Our additional terms depend on two dimensionless parameters, 

2-kGuq 

e= 1 

TK 2 

and a/(/cr), for mass column density a , radius r, epicyclic frequency k, and sound speed 
a. Typically e ~ 0.1, which is small, so our new results are not important modifications to 
previous studies that considered only small Q. However, in regions where Q is large, the 
additional terms lead to residual instabilities that can be important in some situations. 

Numerical and analytical solutions to the modified dispersion relation are found here 
for typical regions in galactic and other disks. These include the main disks of spiral 
galaxies, where the rotation curves are approximately flat (Rubin et al. 1985); the inner 
disks of galaxies, where the rotation curves are approximately solid body; inner solid-body 
gaseous disks that are not self-gravitating (e.g., NGC 2207; Elmegreen et al. 1998), and 
non-self-gravitating Keplerian disks, as might be appropriate for proto-planetary disks or 
galactic nuclear regions surrounding black holes (Nakai et al. 1993). 



2. The General Dispersion Relation 

The dynamical response of an infinitely thin fluid disk to perturbation density 
waves will be studied here, considering various degrees of approximations using algebraic 
expansions in terms of small parameters. The disk response to spiral waves is considered to 
be weak enough for the linearized equations of motion to be valid. The effects of self-gravity, 
pressure, and differential rotation are included. The pressure is assumed to depend only on 
the density; in the formulation, enthalpy is used. In the analysis, perturbation variables 
are assumed to be of the form gi(r, 6,t) = G(r)e l f k ^ dr e t ( ujt ~ me \ where r is the radius, 9 is 
the azimuthal angle, uj is the frequency of oscillation if it is real, and the growth or decay 
rate if it is imaginary, m is the number of arms, k(r) is the radial wavenumber, and G{r) 
is the slowly varying amplitude. The spiral waves have an interarm spacing that is much 
shorter than the radius, that is ( = l/|Ax| <C 1 for total wavenumber k = \Jk 2 + m 2 /r 2 . 
This condition is satisfied either for very short waves or for open spirals with many arms, 
and it allows asymptotic solutions to the density response. The same condition is used to 
express the density as a linear function of the gravitational potential (Bertin & Mark 1979). 

The linearized equations of motion are combined with the continuity equation to relate 
the perturbation enthalpy hi to the perturbation gravitational potential 4>i (Goldreich & 
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Tremaine 1979, Lin & Lau 1979): 



C(h 1 + 4 >1 ) = -Ch 1 , (2) 

where C = d 2 /dr 2 + A d/dr + B and the coefficients are A = — (l/r)d\nA/d\nr , 
5 = - m 2 /r 2 + (2mVL/r 2 Kv)d\Ti(K 2 (1 - v 2 ) /a Q)/dlnr, and C = -k 2 (1 - z/ 2 ) /a 2 ; also 
^4 = k 2 (1 — v 2 ) j (er r), where z/ is the dimensionless frequency, v = {u — rriQ) / k, m is 
the number of arms in the spiral pattern, k is the epicyclic frequency, ctq is the surface 
density of the disk, fl (r) is the angular frequency, and a is the sound speed in the disk. The 
perturbation gravitational potential can be expressed in the form 4>i(r) = $(r)e i ^ fc(r - )a!r ; 
then Poisson's equation is (Bertin & Mark 1979): 



0o 



f(r)<f> 



i- 



(3) 



with the definition 

m 



i 



2irGrX(a, m) 
and the approximation 



1 + iA(a) r—+ B(a) r 2 -— + C (a) (r — 2 
ar c(r dr 



K(a, m) 



T 

a 



m + 1/2 9x 
r(l + ^-^r 2 ), 

( a 2 + ( m + l/2) 2 )" 1/2 , 



(4) 



This expansion for f(r) is correct to third order in a. The terms A, B, and C are defined 
in Bertin & Mark (1979); they are: 

A(a) = X 2 -X x 2 

B(a) = Ki 3 + X 3 - 2AT 1 A: 2 

(7(a) = 9X 1 2 X 2 -QX 1 X 3 + 3X 4 -3X 1 * -3X 2 2 , 



where 

1 d n X(a } m) 



K„ 



n\X(a,m) da r 



The enthalpy, /zi = a 2 cri/o"o, can be expressed in terms of the potential 4>\ using 
equation (|3|) to obtain 

fci = -/(r)&. (5) 
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The expression for /, equation (^), can be expanded in the small parameter ( to get f(r) = 
(k/kj) [1 + i fx ( + f 2 ( 2 + (fz + i U) C 3 + •••]• Here, kj = 2irG(T /a 2 is the two-dimensional 
equivalent of the Jeans wavenumber. The terms /, are real and depend on derivatives of k 
and $. For example, fx = (k/k) —l/2 — r&/$ — (l + rk'/k) (m 2 /2k 2 r 2 ^ . 

If only the first term is kept in the expansion of f(r) and all radial gradients and the 
m-dependence of k is dropped, the Lin & Shu (1966) dispersion relation is obtained: 



(u-mSi) = k 2 - 2nGa Q \k\ + k 2 a 2 . (6) 

In terms of the dimensionless frequency, v = (uo — mVt) / k, dimensionless wavelength, 

V = k crit /\k\ > 0, where k crit = k 2 /27iG(Tq, and Toomre's stability parameter Q = Ka/nGao, 

the Lin-Shu relation is 

7] AT] 2 



3. Tangential forces and the stability parameter J 

3.1. The Bertin-Lin-Lowe-Thurstans dispersion relation 

In the derivation of the Lin-Shu dispersion relation, which is equation (0) above, 
terms of magnitude m/kr are ignored. Thus the dispersion relation is accurate for radial 
oscillations only. When the azimuthal wavenumber m/r is included, the gravitational 
instability is stronger (Lau & Bertin 1978). In the derivation of the corresponding dispersion 
relation, Lau and Bertin made the assumptions that in Poisson's equation the out of phase 
(i.e., imaginary) terms can be ignored and the wavenumber \k\ ~ kj/2. Defining the total 
wavelength to be A m = 2n / \Jk 2 + m 2 /r 2 , Poisson's equation becomes 

- 01 = GcxxXm, (8) 

and equation @ is 



(O~l/Oo) iri pha 



h + 



K 



mVt) 



Tx 



(1 



(9) 



where Tx = —(2mQ/K,r) 2 (d\nQ/d\nr). Note that the last term in the equation above 
contains (1 — u 2 ), which was not present in (C15) in Lau & Bertin (1978) because they 
were considering solutions near corotation [y ~ 1). However, Ti/(1 — u 2 ) can be derived 
from their equations (B6) and (B9), it comes from their second term in equation (B9); in 
fact, Bertin et al. (1989) included it in their dispersion relation. Lau & Bertin (1978) also 
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dropped the fifth term in (C14) when they derived (C15) because it is higher order in 1/kr. 
We do the same for equation fl9|) because this section is about the low order terms as well. 
We include all of these terms in the higher order analysis in the rest of the paper. 

The dispersion relation for spiral waves, which is analogous to equation (|7|), is now 

Q 2 /i = v- --J 1 ,!/^ 2V ( 10 ) 

fj 1 + J /{l — v ) 

where f] = k crit /k and J 2 = Ti/k 2 rit , as defined in Bertin et al. (1989). We call equation fllPf ) 
the Bertin-Lin-Lowe-Thurstan (BLLT) dispersion relation with dimensionless frequency 
v bllt- It describes the response of a differentially rotating disk to spiral perturbations. 
Evidently, the response is stronger than for axisymmetric perturbations by a factor that 
depends on the parameter J. 

Equation ( ]TTJ| ) was studied extensively by Lau & Bertin (1978) and Bertin et al. (1989) 
in the limit when v ~ 0, which is near corotation. In this limit, equation ( |TUD predicts an 
instability when the frequency is purely imaginary, and this occurs when 

J 2 fj 2 ) < 0. (11) 

For ring-like perturbations (m = and J = 0), equation (|TT|) is satisfied when Q < 1; 
that is, equation (|TTD reduces to Toomre's (1964) instability condition, Q < 1, for the 
axisymmetric case. 





It is seen from equations (^) and (10) that when the imaginary terms in the equation 
of motion and Poisson's equation are ignored (Hunter 1983), the dimensionless frequency is 
pure real or pure imaginary according to the values of fj and Q and for small values of J 2 . 
The exclusion of these imaginary terms is justified in the limits \kr\ » 1 and k cr i t r » 1. 
This latter quantity is e -1 , defined by equation ([[]). If the complex terms are included in 
the equation of motion and Poisson's equation, then the frequency solutions are complex 
functions of fj and Q. In that case, the frequency v contains a non- vanishing imaginary part 
in all of the parameter space {t],Q)- This means there is always some instability present, 
consisting of an oscillation plus growth, so Q is not an absolute discriminant of stability for 
small J when higher order terms in e are included. These new instabilities will be discussed 
in detail in sections ^ and |5], but first we consider the low-order BLLT equation in the 
region beyond the Lindblad resonances. 
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3.2. A modification to the BLLT equation beyond the Lindblad Resonances 

In addition to the instability condition given by equation (pi]) , the BLLT dispersion 
relation (Eq. |It]) predicts another instability when the frequency v is complex and has a 
real component with an absolute value larger than 1. 

This is a different regime of position relative to the resonances than considered by 
BLLT. They were concerned mostly with instabilities near corotation, where the waves 
are evanescent. For this reason, they took v ~ 0. In this section, we consider stability 
properties inside the inner Lindblad resonance [y < —1) and outside the outer Lindblad 
resonance (v > 1), using the same order of approximation as in BLLT. These are regions 
that were considered to be damped and radiative, respectively, in the BLLT model. We 
show that the same dispersion relations also allow solutions that grow as they oscillate, i.e., 
with complex frequencies. 

The condition for this second instability may be obtained from the square root part of 
the solution for v 2 in equation ([U]), and is: 

|-i-« 2 r)<o. (i2) 

This condition can be written in the form 

- < % < ~ + (13) 

which is the same as 




f U, IX f <±<5 III, 

ekr < — — < ekr + — (14) 

.-2 TOO \ ' 



a 2 k f As 2 m 2 

— — < ekr H — - 

K 2 k 2 r 2 

if we substitute Qe/2 = a/(nr) and efj = 1/kr, and define J 2 /e 2 = s 2 m 2 , where 
s = 2(— VLrVl) 1 / 2 / k and is of order 1. Equation (|14D is a new condition for instability. When 
this condition is satisfied, the self-gravitating disk is unstable to the growth of spiral waves. 
The right hand side of equation ([14]) contains two terms. The first term depends on the 
self-gravity of the disk and the second depends on shear. When gravity is negligible, there 
is still instability from the second term, coming entirely from pressure, shear, and Coriolis 
forces. We refer to this as an acoustic instability; it has apparently not been considered 
previously in the literature. 

Figure 1 shows the unstable regions for a five-arm spiral (m = 5) in the {k cr it/\k\,Q 2 ) 
plane from the BLLT dispersion relation, equation (|I~0D, considering a self-gravitating disk 
with a flat rotation curve (s 2 = 2); this case is studied in more detail in the next section. 
The growth rate is represented as a gray scale, and the borders of the regions of instability 
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are represented as lines, obtained from the instability conditions. The most unstable region 
is in the bottom left corner of the figure, where the bottom line shows the stability limit for 
the Lin-Shu dispersion relation (m = 0), which is obtained from equation (0). For m and 
J 2 7^ the border of this region of instability shifts to the line given by the BLLT condition 
(Eq. The acoustic instability is bracketed by the two upper lines described by equation 
(|l3l). The lower line corresponds to v 2 = 1 . This occurs at a Lindblad resonance when 
the Doppler-shifted frequency of oscillation, (u — mQ), matches the epicyclic frequency, 
k, which is where the self-gravity of the disk is balanced by the pressure force (the Jeans 
condition) according to equation (|T0|). The upper line corresponds to v 2 = 1 + 2J 2 f) 2 . 

We can investigate the instability conditions (|TTJ) and ( fUj ) further by writing the BLLT 
dispersion relation without self-gravity. This can be done by multiplying equation (PH|) 
by e 2 , and then substituting as above. We then let e — > to turn off gravity. The BLLT 
dispersion relation becomes 



a 2 [k + s 2 m 2 /r 2 



- (2 + a 2 k 2 /K 2 ) v 2 + 1 + — ^ J - = 0. (15) 



We combine the contributions to the dispersion relation from the sound speed and the 
epicyclic frequency by defining an angle 7 = taii" 1 (ak/K). We also define an angle 
p = tt + tan _1 (m/£x) for k < 0; this angle is between tt/2 and tt, giving ship > and 
cosp < 0. The standard definition of a spiral arm pitch angle is tt — p. For e = 0, equation 
(PI) is never satisfied, so the BLLT instability disappears, as recognized by these authors. 
However, the acoustic instability remains, with an instability criterion given by equation 
(|TJ) with e = 0; this is 

a 2s sin p 2sm 

~ < ~t = TTTs 2- 16 

kt kr k^r z + m z 

Another way to write equation ([16]) is to remove the explicit radial dependence; then the 
instability condition becomes 

tan7 = — < 2s sin p. (17) 



The left hand side of the inequality in equation fll7[) is the ratio of the length scale for the 
epicyclic oscillation to the interarm spacing. This ratio has to be less than order unity for 
the instability to develop, which means that there has to be room for epicyclic motions 
within the distance that separates the spiral arms. That is, spiral waves will grow at all 
wavelengths that have enough room for epicyclic motions at the local sound speed. 

When equation ( |T7j ) is satisfied, a non- self- gravitating fluid disk with differential 
rotation will be unstable to spiral perturbations inside the ILR and outside the OLR. For a 
disk with solid body rotation, s — 0, for a flat rotation curve, s = a/2, and for a Keplerian 
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disk, s = v^6, so condition ( |TTD is more easily satisfied, and the growth of instabilities 
is stronger, with greater shear. From equation (|T^), the phase velocity, c p h, and the 
group velocity, c g , of the acoustic waves in the radial direction can be obtained. Define 
z = (2s ship/ tan7) 2 , and w± = (1 ± y/1 — z)/2; then 

Cph = mVl/k ± — sj 1 + w ± tan 2 7, 

w± cos p tan 7 

°g = ±a - r 9 • 

yl — zyl + w±tan 7 

In the unstable regime, z > 1, which implies that c p h, and c g are complex; only the real 
parts should be taken for the physical phase and group velocities. Note that for trailing 
waves, which are the only waves considered here, cosp < 0. This instability, along with 
additional instabilities resulting from higher order terms, will be studied further in section 
|5] for the cases with flat rotation curves and Kepler rotation. 



3.3. Physical insights to the BLLT extension beyond the Lindblad resonances 



The acoustic instability determined by the condition fllTD has a different physical 
origin than the higher-order instability discussed in the next sections. For example, the 
higher-order instability works with or without shear, but the BLLT extension beyond the 



Lindblad resonances requires shear (s 7^ in equation [IT]). We show in section |6| that the 
physical origin of the higher order instability is a geometric growth of incoming wavetrains 
near the nucleus of a galaxy. We do not actually think of this higher-order growth as an 
instability because it is limited in time to the propagation time over the radius. This is 
unlike the acoustic instability discussed in the previous section, which is a true instability. 
The acoustic instability is very similar to the gravity-driven instability of BLLT near 
corotation, i.e., between the Lindblad resonances, but it is pressure-driven instead, and 
beyond the Lindblad resonances. We explain here in physical terms how it works. 

In normal galactic spirals between the Lindblad resonances, and in bars between 
corotation and the ILR, the spiral or bar perturbation grows with time because more and 
more stellar (or fluid particle) orbits lock into phase with the perturbation, and because 
each new aligned orbit reinforces the perturbation, causing greater and greater forcing. This 
works for two reasons: (1) In this radial range, an unperturbed epicycle precesses slower 
than the pattern speed, i.e., the precession speed, Q — n/m, is less than the pattern speed, 
Q p . (2) The inward forcing from the perturbation, gravity in this case, is greatest near the 
apocenter of the epicyclic orbit. For a spiral arm, this apocenter occurs just outside the 
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potential minimum of the arm, and is directed inward because of the arm gravity. For a 
bar, the apocenter is on the bar major axis, and is directed inward because of the gravity 
of the bar. 

Reason (1) implies that in the absence of forcing, a fluid element with its apocenter at 
the crest of one arm will come in and go out again to the next apocenter before it reaches 
the next arm. This is because the precession rate is slow and the apocenter of the epicycle 
twists around in angle more slowly than the spiral pattern. (In other words, the Coriolis 
force (in k) is too large, so the angular velocity perturbation causes too large a radial 
velocity perturbation and the radial oscillation period is short.) However, with gravity, 
the excess inward forcing at the apocenter in the first arm crest gives the orbit an extra 
kick in the radial direction, and this flings the fluid element all the way around to the 
next arm before it has its next apocenter. Moreover, this kick occurs during the part of 
the orbit when the fluid is most susceptible to gaining momentum, i.e., when it is moving 
most slowly and spending the most time (at apocenter). Thus the forced orbit aligns with 
the perturbation, always having its apocenter in the arm crest. The same occurs for a bar: 
the presence of an excess inward bar forcing on the major axis of the bar flings the fluid 
elements around so they have their next apocenter at the other major axis, rather than too 
early. Thus we see how the gravitational force from spirals and bars causes the epicyclic 
motions of individual fluid elements to align with the perturbation and strengthen it. 

Inside the inner Lindblad resonance, the precession speed of an unforced stellar orbit is 
greater than the pattern speed, i.e., Q — k/ui > Q p , so normal spiral or bar gravity kicks the 
stellar orbits the wrong way. That is, the gravity forcing makes an epicycle that already has 
its next apocenter come too late, meet the next arm even later. For the case of the bar, this 
leads to a perpendicular alignment of the orbits, so the point of maximum inward forcing, 
on the bar axis, is at the place in the epicycle, its pericenter, where the fluid element is least 
susceptible to acquire excess momentum, i.e., where it is moving most quickly. A previous 
description of this process was given in Elmegreen (1997). 

Now consider the influence of pressure on these waves. The pressure forcing in a 
spiral is out of phase from the gravity forcing. When the gravity forcing is a maximum 
inward, just outside the spiral potential minimum, the pressure forcing is a maximum in the 
outward direction, because of the pressure gradient from the compression in the arm crest. 
Thus pressure is a stabilizing influence on normal spirals and bars between the Lindblad 
resonances, as is well known. Pressure forcing has to be less than gravity forcing for the 
spiral to grow. This is the usual condition for the dispersion relation, which equates the 
wave oscillation frequency to positive (and therefore stabilizing) contributions from acoustic 
and epicyclic oscillations, plus a negative (and therefore destabilizing) contribution from 
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self-gravity. 

Inside the ILR and outside the OLR, the role of pressure and gravity change. Whereas 
self-gravity opposes the alignment of epicycles beyond the Lindblad resonances, as discussed 
above, pressure is in the right phase to support this alignment of epicycles. The maximum 
outward force from pressure is near the epicycle apocenter both inside and outside the ILR, 
and the existence of this outward force slows down the fluid at its apocenter in both cases 
too. But inside the ILR and outside the OLR, this slow down causes the next apocenter to 
occur in the next arm, rather than after the next arm, which would be the case without the 
pressure forcing. 

The acoustic instability beyond the Lindblad resonances is therefore due to a reversal 
in the role of gravity and pressure as driving agents for spiral density waves on either side of 
the Lindblad resonances. Between the ILR and OLR, gravity changes the orbits in such a 
way that they reinforce an initial perturbation, while pressure opposes this change. Beyond 
the Lindblad resonances, pressure changes the orbits to reinforce the initial perturbation, 
while gravity opposes. When gravity is weak beyond the Lindblad resonances, pressure 
alone is left to drive spiral instabilities. 



The sensitivity of the instability condition fll7|) to shear (s) and pitch angle (tt —p), 
which is the same as the requirement that v 2 < 1 + 2 J 2 r) 2 , makes sense for such pressure 
driven spirals. When the pitch angle is large, the maximum inward pressure force occurs 
closest to the minor axis of the epicycle, and the maximum outward pressure force occurs 
closest to the major axis. This situation leads to the maximum possible forcing from the 
pressure gradients. The shear is important because this is what causes the epicycles to 
precess forward or backward relative to the pattern. Without shear, the precession speed is 
zero, and no amount of pressure forcing can enhance the spiral alignment of orbits. 



4. Higher Order Terms in the General Dispersion Relation 

For the general case with self-gravity, it is possible to solve for the complex frequency 
if we know the basic state of the disk. If the rotation curve, the density distribution, and 
the sound speed distribution in the disk are known, then the dispersion relation in the 
tightwinding approximation can be obtained to second order in e. 

The dispersion relation for v is obtained by turning equation (0) into an algebraic 
expression. This is done by using the definition of the enthalpy and equation (||) to express 
the enthalpy as a function of the potential and then using the asymptotic form of the 
potential, (pi = &e l J k ^ dr . We will consider only trailing spirals {k < 0). Note that 
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v = — {m&i / k + vk I k) for radial derivatives denoted by primes. Equation (Q) can be 
written in the form: 

~r 2 d 2 



— -5- + (Ar — + Br 2 
dr dr 



{h l + <f) l )=5- 1 (l-u 2 )h 



where 5 = a 2 /(K 2 r 2 ). Multiply equation ([18]) by 5 /hi and define D = 5(1 — A), 
D i = 6 fci(<h + hi and ^2 = - m 2 )^! + h). Then 

D 2 + (Ar) D 1 + (Br 2 + m 2 )D + v 2 - 1 = 0. 



(19) 



The terms _D ? - are 
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These terms are used to find numerically the roots of the dispersion relation. They can be 
expanded in the small parameter 1/kr by using the Bertin & Mark expression of Poisson's 
equation (Eqs. ^ and [5|). Their expansion is correct to third order in 1/kr, so our dispersion 
relation is limited to third order in this quantity as well. In terms of Q 2 and rj, and to 
lowest order in e, Dj become: 
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^-v) +te 3 f l 2 f 1 + 
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These equations define the terms d^j\ note that alternate terms are imaginary as is typical 
for WKB approximation methods. Also note that rk' /k = cos 2 prk' /k — sin 2 p, and that 
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rkj /kj = r<7o'/°"o — 2ra /a. We take k to be constant and real. The terms (Ar) and 
(Br 2 + m 2 ) contain contributions of order unity divided by v and [y 2 — 1). To get a 
polynomial expression for u, we calculate the expressions 

v(v 2 — l)Ar = ai v + a 2 v 2 + 03 i/ 3 
z/(z/ 2 - 1) (5r 2 + m 2 ) = b + b 1 u + b 2 u 2 



with 



2f k r cr 

ai = 1 — 



a 2 



2 m Q r Q 



a 3 = 1 + 



ra ' 



00 



2 777, , r OY) ^ ^ 2rs 
&o = C— + 



/2 m fi\ 2 r f2' t9 . 9 9 9 

61 = - (—J IT s J /£ = s m 

2mQ .ran rO , 

62 = — + -0 - - 



Equation ( Jl9| ) is now multiplied by z/(z/ 2 — 1) to obtain a general dispersion relation for 
fluid disks: 

z/ 5 + c 3 z/ 3 + c 2 z/ 2 + c x i/ + c = 0, (20) 

where 

c 3 = -2 + d 20 + ie(d 21 + a 3 dn) + e 2 (d 22 + a 3 d 12 ) + ... 

= c 30 + iec 31 + e 2 c 32 + 
c 2 = iea 2 dn + e 2 (a 2 d 12 + b 2 d 02 ) + ... 

= zec 2 i + e 2 c 2 2 + 
ci = 1 - d 20 + d 02 J 2 + ie(—d 2 i + a x dn + J 2 d 03 ) + ... 

= c U) + i£C 11 + e 2 c 12 + 
c = e 2 Mo2 + ie 3 b d 03 + ... 

= e 2 c 02 + «e 3 c 03 + •••• 
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This dispersion relation includes terms that have been neglected in previous studies. The 
effect of the higher order terms can be followed by the dependence of the coefficients q 
on the small parameter e. In the limit of e — > 0, but with a finite fj and finite Q 2 /r) 2 , the 
general dispersion relation (Eq. |2H| ) becomes the BLLT dispersion relation (Eq. [TOD . 

We investigate the effects of the higher order terms by expressing v as an expansion in 
the parameter e, that is, v = u Q + ev\ + e 2 v 2 + and by solving for the roots of equation 
(|j~9|). Substituting the expansion for v into equation (|T9| ) and setting coefficients of equal 
powers of e to zero, we obtain expressions for the expansion terms vi. The zeroth-order 
root, z/ , satisfies the equation 



z/ 4 - z/o 2 (2 - d 20 ) + 1 - d 20 (l + J 2 fj 2 )] = 0. (21) 



The expression in the squared brackets of equation (|2lD is the BLLT dispersion relation as 
discussed above. The other solution (uq = 0) has no terms of order e; i.e., it is of the form 
v = v 2 e 2 + z/ 3 e 3 + .... 

The first-order term that corresponds to the nonzero solution z/ is 

_ _ { "0 (gu + c 2 i z/q + c 3 i z/p 2 ) ^ 
1 Cio + 3c 30 z/ 2 + 5z/ 4 

for real z^o (i.e., stability in the BLLT equation), this v\ is purely imaginary; for imaginary 
z/ , it is complex. 

The coefficients, Ci 2 , for the next term in the expansion, u 2 , are pure real and, if uq is 
real, then this term is real also, and a factor of e 2 smaller. When vq is real, the growth rate 
to first order in e is attributed to V\. The next contribution to the growth rate will be from 
z/ 3 , which is of order e 2 smaller. 

In summary, we have found in this analysis a general dispersion relation that includes 
the effects of radial variations in the basic parameters of the disk and is accurate to higher 
order in the small parameter e = (/c cr i^) _1 . Furthermore, the effects included in this 
analysis change significantly the criterion for stability of the disk as shown explicitly by the 
models in the next section. 



5. Instability models including the high order terms 

Several models will be studied to illustrate the effects of the higher order terms in 
the dispersion relation and to investigate how different assumptions affect the stability 
of the disk. Four models will be considered: a self-gravitating disk with a flat rotation 
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curve, a self-gravitating disk with solid body rotation, a non-self-gravitating disk with solid 
body rotation, and a non- self-gravitating disk with Keplerian rotation. The amplitude of 
the wave is assumed to be slowly varying so r$'/$ <C 1. This gives an arm/interarm 
contrast that increases with radius beyond one scale length, in agreement with observations 
(Schweitzer 1980; Elmegreen & Elmegreen 1984). 

All disks considered here are assumed to have an exponential mass column density 
profile with a scale length r d and a constant sound speed, a. Then ra '/a = —r/rd and 
r a j a = 0. 

We are considering solutions to the dispersion relation obtained from a local analysis 
where there are gradients in the physical quantities of the equilibrium disk. The local 
analysis is relevant when the growth time for the perturbations is shorter than the time 
needed for the disturbances to travel to the boundaries (e.g., see Lin & Shu 1964; Toomre 
1981). That is usually 10 9 years to the outer boundary and 10 7 years to the center for a 
circumnuclear disk, but in this case the center boundary usually serves as a sink, as waves 
are shocked and energy is dissipated. Therefore we are justified in using a local analysis in 
nuclear disks. For main galaxy disks the growth time of spiral waves is also typically less 
than the propagation time. Bertin et al. (1989) considered a non-local analysis, including 
the effects of gradients and boundary conditions. This leads to the standard modal theory 
of spiral structure. 

Gradients of disk properties, as well as curvature, can lead to spatial variations in the 
amplitude of spiral waves, including singularities. The curvature effects are considered in 
more detail in section || 

General dispersion relations like these can be solved by assuming k real and uj complex, 
or k complex and u real. In the remainder of this section, we consider k real and constant 
and look for solutions with imaginary u. The result will be sinusoidal waves that grow 
exponentially with time, as in the usual stability analyses. 

A third method of analysis is to consider the initial value problem of time- dependent 
growth with shearing sinusoidal perturbations, as in Goldreich & Lynden-Bell (1965) and 
Toomre (1981). When gravity is important, this leads to the swing amplifier theory. 

In the following subsections we will investigate analytically and numerically the 
dispersion relation for disks with different rotational properties. The relevant dispersion 
relation is Eq. (|i~9|). The same dispersion relation with explicit expansion in terms of the 
small parameter e is Eq. (|20| ) for self-gravitating disks. Another dispersion relation is 
derived for non- self-gravitating disks in the indicated subsections. 



-16- 



5.1. Exponential self-gravitating disk with constant rotation velocity 

We first find the roots of Eq. ( |T9"D at two scale lengths for an exponential disk with a 
constant rotation speed. In this case rO /Q = tk/k = — 1, and Q/k = l/V%- The value 
of e = l/(k cr itr) depends on the ratio of the disk to total mass (disk and halo) in the spiral 
region. A value of e ~ 0.11 corresponds to the Solar radius in the Galaxy, using the rotation 
curve model in Schmidt (1983) and a disk mass surface density of 48 M pc -2 (Kuijken & 
Gilmore 1989, 1991). We use a value of e = 0.1. 

There are five roots of the dispersion relation. The root that corresponds to the 
greatest growth is always plotted in the figures here; this is the root with most negative 
imaginary component. 

Figure 2 shows the components of the normalized frequencies v in the (k cr u/\k\ = rj, 
Q 2 ) plane for two values of the azimuthal wavenumber, m =2, and 5, obtained numerically 
using the full dispersion relation, equation (|19|) with coefficients up to third order in the 
small parameter £ = efj. To be clear, we write k crit /\k\ instead of rj in the figures. The top 
figures show the negative of the imaginary component of the frequency, i.e., the growth 
rate normalized to the epicyclic frequency k, with contour values 2*/ 4 for i = -20 to 10. 
The bottom figures show the absolute values of the corresponding real frequencies with the 
same contours. The left figures correspond to m = 2 and the right correspond to m = 5. 
The values of the real and imaginary components are tabulated for some values of Q 2 and 
kcrit/\k\ in table 1; this will facilitate the interpretation of the contours. 

The thick lines in the top plots of figure 2 indicate the loci of points where the 
normalized frequency, vbllt, equals (corotation, lower line), ±1 (inner and outer Lindblad 
resonances, middle line) and ±\/l + 2J 2 7] 2 (upper line) in the BLLT dispersion relation, 
equation (|H]). The BLLT instability condition, equation flTT|), is satisfied below the lower 
thick line. The new acoustic instability condition, equation (114]) , is satisfied between the 
middle and the upper thick lines. 

The figure and table show that the growth rate decreases but remains finite for 
k C rit/\k\ — > 0, and that at k cr u/\k\ = 0, it increases with increasing Q. At intermediate 
values of k crit /\k\, say 0.5, the growth rate is largest for Q < 1 and decreases to a minimum 
at Q 2 ~ 2, but again increases for increasing Q 2 . The growth rate decreases for increasing 
k C rit/\k\ beyond 0.5 for constant values of Q 2 . This pattern is observed for both m values. 
A significant difference between the figures for m = 2 and m = 5 is that for higher m, the 
growth rate is larger over the plotted (k cr i t /\k\, Q 2 ) plane than for low m, and for high 
k C rit/\k\, the growth rate remains relatively large for moderate values of Q 2 above the line 
v bllt — 0. This enhanced growth at high m is because the J-parameter is proportional to 
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m and is contributing to the higher order terms in the dispersion relation. 

Note that there is a kink in the lower right corner (k cr u/\k\ « 1.6) of the m =2 contour 
plot for the real component of the root. This occurs because in adjacent regions to the kink 
different real components have the most negative imaginary component. 

Figure 2 and table 1 also indicate that the greatest growth occurs for small values of 
Q 2 , just as predicted using the BLLT dispersion relation (cf. Sect. pT2| ). Moreover, they 
indicate that the disk is unstable to form spirals for a wide range of Q and m, although 
the growth rate is low, of order e, when Q is large. This implies there is still a spiral 
instability at low gravity. For most bright galaxies, however, the region where the rotation 
curve is flat is also the region where Q is relatively small, so these high Q solutions are not 
important. They could be important in early type galaxies (Caldwell et al. 1992) or low 
surface brightness galaxies (van der Hulst et al. 1993) where Q is high in the main disk. 



5.2. Exponential self-gravitating disk with solid body rotation 

The inner parts of galaxies and small galaxies typically have rotation curves that are 
approximately solid body. This is the result of a strong bulge with a nearly uniform central 
density in some spiral galaxies, and a relatively dense dark matter halo in dwarf galaxies. 
Inner galaxy disks (Elmegreen et al. 1998) and dwarfs (Hunter et al. 1998) may also be 
weakly self-gravitating for some time (e.g., between accretion events and starbursts), and 
so the high-Q cases studied here may have applications there. Furthermore, inner disks and 
dwarfs have short rotation times, so the actual growth factor of a spiral instability can be 
large even if the normalized growth rate is small. 

For solid body rotation, rfi'/fi = tk/k = 0, and Q/k = 1/2. We assume a value for 
e = 0.1 as in the previous section. In this case the term A(r) does not depend on v and 
B(r) has a 1/u dependence. The dispersion relation then becomes cubic in v: 

+ {-l + D 2 + azD 1 ) u+ b 2 D Q = 0, (23) 

where the terms Di, a^, and b 2 were defined in the previous section. The roots can be 
expressed as an expansion in e, writing v = u + V\ e + v 2 e 2 + .... The zero order term is the 
Lin-Shu dispersion relation, equation (H), with rj replaced by fj. The first order term is 



. cosp 
V\ e = — i e 



Q L ra . 2 \ 1 ra sin p 
— 1 2 sin p 



(24) 



2u 

In the region where v is real, the growth rate is dominated by the first order term. In the 
region where \uq\ is of order 1, V\ ~ —i and the growth rate is of order iv\ e ~ e. For an 
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exponential disk at two scale lengths = —2, so 

(3-2 sin 2 j9) + 



Vi e = — 1 e 



CO 

cosp 



2u 



Q 2 / 9 \ sin 2 p + 3 

^(3-2sin 2 p)+ P - . (2:)) 



Figure 3 and table 2 show real and imaginary components of the normalized roots of 
the full dispersion relation (|23| ) for the rising rotation curve model at r = 2 r^. Again we 
display only the root that corresponds to the fastest growth. We can see from the left-hand 
regions in the (k crit /\k\, Q 2 ) plot, where the absolute values of the real components are 
large, that the growth rates become small for small k cr u/\k\. The opposite occurs for small 
values of the real component, which are in the lower region of the plot. Where the real 
component is of order 1, in the center of the plot, the growth rate is of order e. 

The detailed behavior of the growth rate in this case can be followed from the 
approximate analytical solution written above as equation (p4|). For example, equation 
(p4j) gives the same growth rate as the full solution in table 2 for 7] = k crit /\k\ = 0.2 for 
both m = 2 and 5, because the approximate equation is relatively accurate for low fj. 
Equation (23) gives slightly different rates than table 2 for rj = 0.6; at m = 2 and Q 2 = 2, 



5, and 10, equation (p^) has growth rates of 0.232, 0.228 and 0.275 while table 2 has more 
precise growth rates of 0.227, 0.225 and 0.273. The rates given by equation (^) differ more 
significantly from those calculated by equation fl2"3]) when 77 > 0.6. 



One can observe from table 2 that the real component corresponding to the greatest 
growth rate in the (k cr i t /\k\, Q 2 ) plane is always negative, i.e., it corresponds to the Lin-Shu 
and BLLT solutions inside of corotation in the disk. 

Figures 2 and 3 show that there is a similarity between the growth rates for the flat 
and solid body rotation curve models. Both figures display a saddle shape for the growth 
rate contours; the greatest growth occurs as Q 2 < 1, and for k crit /\k\ ~ 0.5, the growth rate 
first decreases and then increases with increasing Q 2 . The main difference between the two 
models occurs for large numbers of spiral arms, where the growth rate is smaller at m = 5 
than m = 2 for the solid body case, and larger at m = 5 than m = 2 in the flat rotation 
curve case. This is because for solid body rotation, J = 0, so the absence of differential 
rotation reduces the growth rate of waves at any m^0. The zero order BLLT instability 
condition (Eq. (|TI|)) is reduced to the Lin-Shu instability condition (Eq. (|7j)), and the 
acoustic instability disappears as the upper unstable region collapses around v 2 « 1. In 
addition, the contributions of J 2 to the higher order terms 2^ are also absent so the growth 
rate is less than for J 2 > 0. 

The solutions shown for all the self-gravitating models indicate that disks are weakly 
unstable to spiral waves when e = l/(k crit r) > 0, even in the limit of weak self-gravity. 
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This is the first time spiral disk instabilities have been found at large Q in the absence 
of magnetic fields. We pursue this result further in the next section, which considers the 
growth of waves in the absence of self-gravity, that is, when e = 0. 



5.3. Exponential disk with solid body rotation and no self-gravity 

This section and the next consider fluid disks without self-gravity as an idealization 
of the high Q cases found to be unstable in the previous sections. To be consistent with 
the radial dependence of the enthalpy amplitude, H(r), used before, which was defined in 
terms of a slowly varying potential amplitude we now assume H(r) oc —f(r), where / 
was given in the discussion following equation (^). 

From equation ([24]) we can see that the growth rate of the instability depends on both 
the self-gravity of the disk and the radial derivative of the background surface density. The 
normalized growth rate is, to first order, v\ e, from the previous discussion. The first term 
in equation ( [24] ) is proportional to Q 2 cospe/(4fi) = ka 2 /K 2 r, which is independent of the 
self-gravity of the disk. It depends primarily on the disk curvature, i.e., on the ratio of 
the square of the semimajor axis of an epicycle caused by random motions (o/k), to the 
product of the wave scale (& _1 ) and the disk radius (r). The second term is proportional to 
e oc massd / masstotai , which comes from the self-gravity of the disk. If the disk self-gravity is 
neglected, e = and the second term is zero, but there is still growth from the first term, 
depending on orbital curvature. When e = 0, the expansion has to be made in terms of the 
small parameter £ = 1/kr. Then we get: 

u = ±\jl + {ak) 2 /K^ 

k Qj^ f 7* \ 

u ^ = l ^~2 2sin 2 p 1 . (26) 

2n z rv V r d J 

When the expression inside the parenthesis of equation ( ^6|) is zero for some particular 
pitch angle 7r — p, there is no growth at that radius, but there is growth at adjacent radii. 

The numerical solutions to equation (^) when gravity is neglected are shown for 
r = 2r d in figure 4 and table 3, using normalized axes (a/ftr) 2 instead of Q 2 and l/|Ar| 
instead of k crit /\k\. To compare the growth rates with the previous models, recall that the 
value of the vertical axis in figure 4 is obtained by multiplying the value of the vertical 
axis in our previous figures by e 2 = 0.01, and the value of the horizontal axis in figure 4 
is obtained by multiplying the previous value of the horizontal axis by e=0.1. This means 
that the growth rates in figure 4 are analogous to those in the upper right part of figure 3. 
Figure 4 and table 3 indicate that the growth rate remains finite, proportional to (a/«x), as 
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\kr\ — > oo. We infer from this behavior that the instability is acoustic in nature, similar to 
that described in section |3]2|, but in the absence of self-gravity and shear. It is driven by 
curvature and pressure gradients in the disk (cf. section IH). 



5.4. Exponential disk with Keplerian rotation and no self-gravity 

Accretion disks around black holes (Nakai et al. 1993) and protostars have negligible 
self-gravity and may have Keplerian rotation. In this case rn /k = rO /Q, = —3/2, and 
equation (0) becomes a fifth order polynomial in v, as for a flat rotation curve. As in the 
previous section, an acoustic instability is still present even in the absence of self-gravity. 
The dispersion relation for this acoustic instability is now obtained from equation (fl9|) with 
the modifications 

D = 5 

D 2 = 5 (-k 2 r 2 + 2ikr^r + T j 1 



/ / J 

Because there is no self-gravity, this fifth order dispersion relation has to be expanded to 
successive orders in £ = 1/kr instead of e, giving v = uq + i/\C, + v^C^ .... The zero-order term 



in this expansion is the modified BLLT dispersion relation, equation (0). Recall that the 
geometric term for a Keplerian disk is s = y6. The zero order term becomes complex when 
the instability condition, equation (|17|), is satisfied. When equation ([H]) is not satisfied, the 
growth rate is dominated by the first order term 

a>k (l + r/r d - 2 sin 2 p) (u 2 - 1) - 3 (1 + mu ) 

viC = 7 x • (27) 

2« 2 r v (2v 2 - 2-a 2 k 2 /K 2 ) 

Figure 5 and table 4 display numerical solutions to the fifth order polynomial, equation 
(|j~9l), for the dispersion relation in this Keplerian model using the modified expressions Di 
when self-gravity is neglected. The real and imaginary components of the root with the 
largest growth rate are plotted using the same axes as in the previous section, (l/|A;r|, 
a 2 /n 2 r 2 ). The critical curve for stability, equation is plotted in the top figures as a 
thick line. To the left of the critical curve, equation (|T6|) is not satisfied and the disk is 
stable against acoustic instabilities to lowest order (higher order instabilities remain). To 
the right of the critical curve, equation (Tl6| ) is satisfied and the acoustic instability to all 
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orders dominates the growth of perturbations. The growth rate is larger than that in the 
case of solid body rotation without self-gravity because shear stimulates growth. There are 
discontinuities in the contours for the real component, with kinks at the same locations 
in the contours of the imaginary growth rate near the critical curve. To the left of these 
discontinuities, the real part of v is negative, corresponding to radii inside the ILR; to the 
right, the real part is positive, corresponding to radii outside the OLR. 



Equation (27]), which is the first order approximation to the growth rate, matches the 
full numerical solutions in figure 5 and table 4 to two significant digits for kr > 5 and 
(a 2 /K 2 r 2 ) < 0.1. 



6. Physical Insights to the Curvature Terms 

We have just shown that differential rotation, curvature, and radial gradients in the 
basic properties of a fluid disk affect the propagation and growth of spiral disturbances. 
Here we simplify the problem by including only the effects of orbital curvature. 

The curvature terms can be illustrated by considering an ideal disk with solid-body 
rotation, constant surface density, and negligible self-gravity. Such disks may be appropriate 
for the central regions of quiescent galaxies, such as NGC 2207 (Elmegreen et al. 1998). 
The governing equation (0) for such a disk, assuming constant surface density, becomes 



d 2 hi 1 dhi I k 2 2 
dr 2 r dr \a 2 




0. (28) 



This equation was derived with the center of the coordinate system at the center of 
rotation. It is the well-known Bessel equation, and can have mathematical singularities 
at r = 0. Often in wave equations, these mathematical singularities can be transformed 
away by a change in the coordinate system, adopting, for example, a rectilinear coordinate 
system instead of cylindrical. However, in the case of a galaxy, the singularities cannot be 
transformed away by a different coordinate system: rotation and galactic gravity define the 
coordinate system. 

In the galactic Bessel equation, the time derivative in the equation of motion appears 
as the term z/ 2 , as it did in the previous sections. The spatial variation in the azimuthal 
direction is also assumed to be the same as before, e~ tme , but in the radial direction it is 
written explicitly. We may look for the behavior of spirals by assuming radial solutions 
of the form hi oc e lkr . These solutions are trailing spirals when k < 0. They always 
contain pieces of waves that can come close to the origin, depending on their direction of 
propagation, so they can force out the singularities in the Bessel equation. The pure-ring 
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case with m = may also approach the origin and increase in amplitude. In this case the 
increase is analogous to laboratory sono luminescence, in which sound waves converge to the 
center of an air bubble in a liquid and increase in amplitude until they shock and emit light 
(Kondic, Gersten, & Yuan 1995). 

In the case of spiral solutions, the radial derivatives in the Bessel equation are replaced 
by ik and the frequency v may be solved to give 



v 1 



— ) 2 (r 2 k 2 + m 2 ) l 

KT 



irk. (29) 



This frequency is necessarily complex because of the first derivative term in equation 
Because of this term, the general solutions are growing or decaying oscillations with spiral 
shapes having m arms. It will become apparent shortly that the incoming waves are 
growing, and the outgoing waves are decaying, as expected from the nuclear singularity. 



For ur >> Vi, we recover the same result as equation flzq) in the limit — > oo, 
assuming constant enthalpy amplitude: 



v R = ±i 



-2.)>*» + m») + x"\ ■* = -(-)*-■ (30) 

KT J v ' 2 \KtJ Ur 



1 / a \ 2 rk 



This gives the growth rate of a wave with azimuthal wavenumber m and radial wavenumber 
k. Note that \ur\ > 1 in all cases here, which means the waves are only inside the inner 
Lindblad resonance or outside the outer Lindblad resonance. This was the case also for the 



new instability solutions discussed in section |3.2| , which re-considered the BLLT equations 
in this new radial limit. 



The nature of the growth implied by vj in equation ( |30|) should be discussed 
more. Recall that the assumed time behavior of the wave in an inertial frame is 
e i{kr+ujt-me) _ e i(kr+uKt) f Qr y _ ^ _ m ^ j K anc j q = Qf following the rotation. Also note 

that we have written v = vr + ii/j. Thus we have a time behavior e %v R Kt e - v i Kt _ When uj is 
negative, the wave grows in time. This occurs for trailing waves only when Vr is the negative 
one of the two solutions given above, because vj cx —k/vR, and k < for trailing waves. 
Moreover, the negative vr solution is an incoming trailing wave, because the wave-like part 
of the solution, e l ( kr + u R Kt ) ; has constant phase for decreasing r with increasing time when 
Vr < 0, i.e., r = — (vr/Is) nt = —\uR/k\Kt when both Vr and k are less than zero. As a 
result, the galactic Bessel equation has trailing spiral wave solutions that grow in time as 
they propagate toward the center of the galaxy. 

These solutions are not instabilities in the usual sense, because t cannot be allowed 
to go to infinity. The waves reach the center in finite time, i.e., in the time t ~ r/a. In 
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this sense, the growing solutions are like those in the galactic swing amplifier (Goldreich & 
Lynden Bell 1965; Julian & Toomre 1966), in which spiral waves grow in the shearing part 
of a disk for a finite time (At ~ 2/ A for Oort constant A). The instabilities in a galactic 
nucleus are also not stationary waves that grow in amplitude without any change in shape. 
This is because the spiral solution is always undefined at the nucleus and can never be 
considered present at all radii. The waves are only pieces of spirals, moving inward or 
outward with a growth or decay in time following the wave crest, respectively. Thus the 
growth is also unlike the growth of infinite plane waves in a sheet, as might be the case for 
the Kelvin-Helmholtz instability, for example. Spiral wave growth in galactic nuclei involves 
inward propagation of finite wave trains. 

The dispersion relation may also be regarded as an equation for k, in which v is 
held as a real variable. Then equation Q2"8| ) has normal Bessel function solutions J m {kBf) 
and Y m (kBr) for 



K 



2 



— > 0; 



ks is the radial wavenumber. When fc^r is large, J m and Y m behave like sines and cosines, 
which may be combined as outgoing or incoming waves with exp(iut). When k B r is small, 
the Y m solutions grow algebraically with decreasing r. The growth arises directly from 
the curvature terms, namely, the first derivative term and the m 2 /r 2 term. These are the 
same terms that led to imaginary v in equations (|29]) and (|30D . The m 2 /r 2 term actually 
defines a region within which the waves start to grow out of bounds, i.e., when r < m/ks 
for 7717^0, the Y m (kBr) solution begins to increase. In terms of the growth discussed for 
the time-dependent case, this is the radius at which an incoming wave has only one more 
epicycle in time before it reaches the nucleus propagating at the sound speed. 

What happens to a trailing spiral wave in a real galaxy when it enters the rfcg/m < 1 
regime? We expect that the amplitude will begin to increase geometrically until nonlinear 
and dissipative effects come into play. This means that the waves will break in the form 
of shocks shortly after they enter the inner region. The condition rks/m < 1 implies that 
the radius for this wave shocking increases with azimuthal wave number m. This explains 
for the case of NGC 2207 (Elmegreen et al. 1998) why the multiple-arm features are only 
observed in the outer part of the nuclear disk, while the m = 1 and ring-like feature is 
close to the center. That is, the multiple arms (high m) become non-linear and damp 
out before they reach the inner radii, leaving only the low-m arms near the center. Other 
spirals that may travel outward in NGC 2207 are probably too weak to be seen because 
their amplitudes decrease as they propagate. 

Galactic nuclear spiral waves also propagate in the azimuthal direction with angular 
speed w/m as long as v 2 — 1 > 0. This angular speed implies that waves with different m 
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will interact, forming complex structures. The waves are also dispersive, with dispersion 
relation 

uj-mVL ( a 2 k B 2 \ 1/2 

-^ = ± ( 1 + ^j ■ 

They form wave packets that propagate with group velocity 

c g = ±asin7 B , 

with = tan -1 aks/n. Undoubtedly the waves will interact because of these various phase 
and group speeds. They will also get sheared by differential rotation in reality to form 
complex spiral structures. When v 2 — 1 < 0, the entire disk is evanescent; then we should 
not see any waves. 

So far we have ignored the exponential density distribution of the disk. If it is taken 
into consideration, equation (123) will change to 



d 2 h x 1 / r\dhi 



dr 2 r \ TdJ dr 



k 2 / o \ m? m 



/ 2 j\ 
n 2 V / r 2 



VTTd 



h = 0, (31) 



where is the scale length of the exponential disk. The additional factor in the first 
derivative term will modify the behavior of the Bessel functions when r>r^, and the 
additional term in the last parenthesis will complicate the wave behavior. But the 
qualitative nature of the Bessel function solutions does not change. 



7. Summary 

We have obtained dispersion relations for spiral waves with multiple arms, considering 
curvature and gradient terms that were ignored in previous derivations. These dispersion 
relations suggest the presence of several new instabilities. Four specific cases were studied, 
flat and rising rotation curves with self-gravity, rising rotation curves without self-gravity, 
and Kepler rotation curves without self-gravity. These cases seem to have applications in 
various regions of galaxies and accretion disks. 

When self-gravity is present, instability at lowest order in the parameter e (cf. Eq. 
[I]) is driven by both shear and self-gravity. Then there are two independent instability 
conditions, either of which can cause spiral waves. These are equations (|TT|) and flU]). The 
first of these comes from Bertin et al. (1989), and contains the Toomre (1964) instability 
condition, Q < 1, as a special case for ring-like perturbations (m = 0, which gives J = 0). 
This first instability is the spiral instability that is commonly discussed in the literature as a 
source of multiple arm and grand design spiral structure in galaxy and protoplanetary disks. 
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The second of these conditions arises outside the Lindblad resonances from a combination 
of parameters different than the first when Q 2 > Afj (cf. Eq. |13]). When self-gravity 
is not present, this second case is still unstable as a result of pressure and differential 
rotation alone, as determined by the smallness of the parameter a/(nr) (cf. Eq. [TBI). This 



pressure-rotation instability is apparently new, and we call it an acoustic instability. A 



physical explanation for it was given in section 3.3 



We also found additional instabilities coming from higher order terms in an expansion 



of the dispersion relation (|19|) around the small parameter e. These additional instabilities 
are present even when the BLLT and Toomre instability conditions are not satisfied, i.e., 
when the low order terms give stability. The source of these residual instabilities is a 
combination of orbital curvature [terms of order l/(Ar)], self-gravity (terms of order e), 
and various disk gradients (ra'/a, ra'/a, etc.), including shear (the J or s terms). Growth 
rates for these residual instabilities were given to all orders in e for flat and rising rotation 
curves by figures 2 and 3 and tables 1 and 2, and they were given to first order in e by 
equation (El) for solid body rotation. The residual instability that arises from self-gravity 



and orbital curvature (through e), discussed in sections 54. and 5T2 , will be called a 
gravitational- curvature instability. The residual instability that arises from a combination 
of pressure and orbital curvature [through a/(/cr)], discussed in sections |573| and pT4| , will be 



called an acoustic-curvature instability, because it operates even without self-gravity. 

These three new instabilities should be important for fluid disks with negligible or 
weak self-gravity, including proto-planetary disks, gaseous disks around black holes, some 
galactic nuclear disks, low surface brightness galaxy disks, and some dwarf galaxies. In these 
cases, zero-order acoustic and higher-order acoustic-curvature and gravitational- curvature 
instabilities can lead to the growth of spiral or other structures in about an orbital time. 
They are most important in the region close to the center where the orbital time is small. 

Non-linear effects arising from these waves may ultimately lead to visible dust lanes 
(Elmegreen et al. 1998) and associated gaseous shocks (Roberts 1969) in even the most 
weakly self-gravitating disks, with the possibility of heightened self-gravity and star 
formation in some of the compressed regions (e.g., Elmegreen 1994). Non-linear effects 
might also promote accretion flows (e.g., Larson 1990). Indeed, the ubiquity of acoustic 
waves in disks implies that galactic nuclear accretion should occur in a wide variety of 
environments with or without shear, self-gravity, or magnetic fields. 
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Table 1. Flat rotation curve at r = 2r d 



v \k\ - 


0.2 


0.6 


1.0 


1.4 


1.8 


0.2 


0.6 


1.0 


1.4 


1.8 




m = 2 


growth rate 








m = 2 


frequency 






10 


0.280 


0.332 


0.340 


0.325 


0.302 


-7.660 


-2.547 


-1.644 


-1.318 


-1.164 


5 


0.218 


0.272 


0.264 


0.236 


0.214 


-5.229 


-1.709 


-1.156 


-0.973 


-0.892 


2 


0.175 


0.236 


0.213 


0.184 


0.168 


-2.924 


-0.864 


-0.673 


-0.658 


-0.664 


1 


0.188 


0.383 


0.221 


0.163 


0.147 


-1.513 


-0.297 


-0.402 


-0.510 


-0.568 


0.1 


1.867 


0.875 
to = 5 


0.410 
growth rate 


0.179 


0.130 


-0.042 


-0.054 
TO = 5 


-0.116 
frequency 


-0.301 


0.457 


10 


0.319 


0.421 


0.453 


0.446 


0.427 


-7.695 


-2.653 


-1.797 


-1.493 


-1.348 


5 


0.258 


0.364 


0.358 


0.309 


0.262 


-5.255 


-1.787 


-1.256 


-1.058 


-0.942 


2 


0.215 


0.275 


0.322 


0.358 


0.385 


-2.941 


-0.853 


-0.529 


-0.420 


-0.372 


1 


0.224 


0.581 


0.559 


0.544 


0.540 


-1.531 


-0.238 


-0.229 


-0.241 


-0.249 


0.1 


1.914 


1.057 


0.810 


0.714 


0.673 


-0.044 


-0.073 


-0.116 


-0.156 


-0.183 



Table 2. Solid body rotation at r = 2 r d 



Q 2 , 


W = 0.2 


0.6 


1.0 


1.4 


1.8 


0.2 


0.6 


1.0 


1.4 


1.8 




rn = 2 


growth rate 








m = 2 


frequency 






10 


0.254 


0.273 


0.268 


0.250 


0.229 


-7.658 


-2.540 


-1.634 


-1.311 


-1.163 


5 


0.193 


0.225 


0.223 


0.206 


0.188 


-5.228 


-1.703 


-1.155 


-0.992 


-0.930 


2 


0.154 


0.227 


0.211 


0.183 


0.164 


-2.922 


-0.881 


-0.723 


-0.727 


-0.748 


1 


0.173 


0.376 


0.244 


0.184 


0.158 


-1.512 


-0.389 


-0.508 


-0.610 


-0.673 


0.1 


1.838 


0.789 


0.361 


0.206 


0.158 


-0.052 


-0.121 


-0.306 


-0.483 


-0.594 






m = 5 


growth rate 








m = 5 


frequency 






10 


0.252 


0.251 


0.222 


0.187 


0.157 


-7.692 


-2.636 


-1.770 


-1.464 


-1.323 


5 


0.192 


0.208 


0.190 


0.166 


0.144 


-5.252 


-1.766 


-1.227 


-1.056 


-0.983 


2 


0.152 


0.215 


0.207 


0.186 


0.169 


-2.938 


-0.901 


-0.709 


-0.680 


-0.674 


1 


0.171 


0.406 


0.311 


0.257 


0.226 


-1.523 


-0.382 


-0.451 


-0.507 


-0.537 


0.1 


1.843 


0.829 


0.487 


0.372 


0.320 


-0.053 


-0.135 


-0.292 


-0.392 


-0.443 
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Table 3. Solid body rotation with no gravity at r = 2 r d 



(—) 2 


nr-r = 0.2 

fcr 


0.6 


1.0 


1.4 


1.8 


0.2 


0.6 


1.0 


1.4 


1.8 




m = 2 


growth rate 








m = 2 


frequency 






1.0 


1.205 


0.482 


0.376 


0.433 


0.462 


-5.584 


-2.816 


0.777 


0.917 


0.995 


0.5 


0.833 


0.311 


0.308 


0.355 


0.377 


-4.024 


-2.160 


0.671 


0.769 


0.820 


0.2 


0.499 


0.168 


0.139 


0.137 


0.131 


-2.673 


-1.608 


0.493 


0.559 


0.598 


0.1 


0.328 


0.101 


0.051 


0.032 


0.022 


-2.026 


-1.355 


-1.294 


-1.279 


-1.274 


0.01 


0.058 


0.014 


0.007 


0.004 


0.003 


-1.148 


-1.045 


-1.037 


-1.035 


-1.034 






m = 5 


growth rate 








m = 5 


frequency 






1.0 


0.684 


0.179 


0.097 


0.067 


0.051 


-7.173 


-5.420 


-5.271 


-5.230 


-5.214 


0.5 


0.474 


0.122 


0.066 


0.045 


0.035 


-5.147 


-3.941 


-3.840 


-3.812 


-3.801 


0.2 


0.287 


0.071 


0.039 


0.026 


0.020 


-3.375 


-2.657 


-2.597 


-2.581 


-2.575 


0.1 


0.192 


0.046 


0.025 


0.017 


0.013 


-2.509 


-2.036 


-1.997 


-1.987 


-1.983 


0.01 


0.039 


0.008 


0.004 


0.003 


0.002 


-1.252 


-1.162 


-1.156 


-1.154 


-1.153 



Table 4. Keplerian rotation with no gravity at r = 2 



(—) 2 


ttt-t = 0.2 

fcr 


0.6 


1.0 


1.4 


1.8 


0.2 


0.6 


1.0 


1.4 


1.8 




m = 2 


growth rate 








m = 2 


frequency 






1.0 


1.606 


1.088 


1.266 


1.365 


1.414 


-5.736 


-2.960 


2.170 


2.081 


2.040 


0.5 


1.220 


0.914 


1.076 


1.139 


1.171 


-4.161 


2.010 


1.818 


1.748 


1.715 


0.2 


0.861 


0.739 


0.821 


0.856 


0.875 


-2.787 


1.621 


1.489 


1.439 


1.415 


0.1 


0.665 


0.597 


0.645 


0.667 


0.679 


-2.121 


1.419 


1.320 


1.281 


1.262 


0.01 


0.249 


0.236 


0.242 


0.246 


0.248 


-1.197 


1.102 


1.068 


1.055 


1.048 






m = 5 


growth rate 








m = 5 


frequency 






1.0 


1.413 


0.961 


1.097 


1.140 


1.160 


-7.262 


4.327 


4.007 


3.885 


3.821 


0.5 


1.214 


1.266 


1.355 


1.382 


1.395 


-5.220 


3.253 


3.076 


3.008 


2.972 


0.2 


1.041 


1.231 


1.276 


1.291 


1.298 


-3.450 


2.394 


2.291 


2.250 


2.228 


0.1 


0.932 


1.089 


1.117 


1.126 


1.130 


-2.607 


1.960 


1.888 


1.859 


1.843 


0.01 


0.506 


0.539 


0.544 


0.546 


0.547 


1.370 


1.229 


1.205 


1.195 


1.190 
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Fig. 1. — Regions of instability in a self-gravitating disk with constant rotation for 5 arm 
spirals (m = 5). The lowest (thin) line borders the instability condition obtained from the 
Lin-Shu dispersion relation (Eq. ([/])) for m and J 2 = 0. The three upper lines bracket 
regions of zero-order instabilities obtained in the Lau-Bertin dispersion relation. The border 
of the Lin-Shu region of instability shifts to the lower thick line for nonzero J 2 . This line 
is from Eq. ( |Tl~| ) for the usual Lau-Bertin stability condition; the two upper thick lines are 
the boundaries enclosing the acoustic instability according to Eq. (|I3]). The shading of the 
unstable regions gives an indication of the growth rate. The most unstable region is in the 
bottom left corner of the figure. 



Fig. 2. — Contours showing the maximum growth rates (top) and corresponding frequencies 
(bottom) for instabilities to all orders of the small parameter e in a self-gravitating disk 
with a constant rotation velocity and an exponential density profile, evaluated at two scale 
lengths. Solutions for two arm spirals (m = 2) are on the left, and for 5 arm spirals (m = 5) 
are on the right. The thick lines are obtained as in Fig. 1., they border the regions of 
zero-order instability for the Lau-Bertin dispersion relation. The bottom thick line is from 
Eq. ( |TTD for the usual Lau-Bertin stability condition; the middle and upper lines are the 
boundaries enclosing the acoustic instability according to Eq. fllcf) . 



Fig. 3. — Same as Fig. 2 for a self-gravitating, exponential disk at two scale lengths, but now 
with solid body rotation. The thick lines on the top figures border the regions of instability 
for the zero-order Lau-Bertin dispersion relation in the case with no shear (J = in Eq. 



Fig. 4. — Same as Fig. 3, but to all orders of the small parameter l/|£x| in the absence of 
self-gravity. All of the unstable growth exhibited in these solutions is from high order terms. 



Fig. 5. — Growth rates and frequencies for a non-self-gravitating disk, as in Fig. 4, but with 
Keplerian rotation. The thick lines on the top figures border the regions of instability for 
the zero-order, Lau-Bertin dispersion relation in a Keplerian disk without self-gravity (Eq. 
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